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Abstract 

A box-ball system with more than one kind of balls is obtained by the gener- 
alized periodic discrete Toda equation (pd Toda eq.). We present an algebraic 
geometric study of the periodic Toda equation. The time evolution of pd Toda 
eq. is linearized on the Picard group of an algebraic variety, and theta function 
solutions are obtained. 



1 Preface 

A cellular automaton (CA) is a discrete dynamical system in which the depen- 
dent variables take on a finite set of discrete values. Although CAs evolve in 
time by simple evolution rules, they often show very complicated behaviour 1 . 

The box-ball system (BBS) is an important CA in which finitely many 
balls move in an array of boxes under a certain evolution rule [2j |3]. This 
discrete dynamical system is obtained from a discrete analogue of the Toda 
equation through the limiting procedure "ultradiscretization" , and displays the 
behavioural characteristics of nonlinear integrable equations [H[S]. In fact, the 
BBS has soliton-like solutions and a large number of conserved quantities [6l [7] . 
Moreover, the solution of the initial value problem of the periodic box-ball sys- 
tem (pBBS) with one kind of ball has been obtained by ultradiscretizing the 
solution of the periodic discrete Toda equation (pd Toda). 

In 1999, Tokihiro, Nagai and Satsuma pointed out that the pBBS with M 
(M > 1) kind of balls (and with capacity one) is obtained by ultradiscretizing 
the generalized periodic discrete Toda equation [8]. Tokihiro and the author 
performed the ultradiscretization of theta function solutions for M — 1 of the 
pd Toda [9]. It is to be expected that the solution of the initial value problem 
of the pBBS with M kinds of balls can be obtained from the solution of the 
generalized Toda equation, as was the case for M = 1. 
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In this paper, we study the generalized pd Toda (M > 1): 
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with periodicity condition I* n+N = /* , KJ' +Ar = V*, where N,n,t £ N. 

The time evolution of this system is linearized on the Picard group of an 
algebraic variety. For example, the time evolution of pd Toda with M = 1 is 
linearized on Jacobi varieties of hyperelliptic curves [TU]. For general M, the 
algebraic curves which appear in the linearization are more complicated than 
hyperelliptic curves. 

These algebraic curves were analyzed by P. van Moerbeke and D. Mum- 
ford in 1978 [11]. In their work, they clarify the algebro-geometric properties 
of the curves which are defined by so-called regular periodic operator matri- 
ces. Although our matrices are not always regular operators, their results are 
essentially also applicable to our case. 

The special feature of the discrete system (|1.HI1.4[) is that we can explicitly 
determine the action of the unit time evolution ( i-> i + 1 on the Picard group 
( proposition 1 2 . 1 6[ ) . Using this result, we extended the theta function expression 
of the solution of pd Toda to the M > 1 case. 

The paper is organized as follows. In section[2j we introduce a spectral curve 
associated to the discrete system (|1.HI1.4[) and discuss its algebro-geometric 
properties. Our aim in this section is to determine the actions of the index shift 
n i— ► n+1 and the time shift t^f+lona Picard group Pic d (C) of this spectral 
curve C. In section[31 we give the theta function expression ftheorem !3.2p which 
is the extension of the formula obtained by Kimijima and Tokihiro 10J . 



2 Spectral curve associated with pd Toda 
2.1 The nature of the spectral curve 

The pd Toda equation (jTTTJ) , (|1.2[) is equivalent to the following matrix form: 



L t+1 (y)R t+M (y) = R t (y)L t (y), 
where L t (y) and Rt(y) are given by 

vi 1 

L t (y)= . . My) 

1 / 
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v N-l 



(2.1) 
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and y is a complex variable. Let us introduce a new matrix X t (y) denned by 

X t (y) :=L t (y)R t+M ^(y)---R t+1 (y)R t (y). (2.2) 

From (j2~Tj) and (j!T2")) . we obtain 

X t+1 (y)Rt(y) = Rt(y)X t (y), (2.3) 

which implies that the eigenvalues of Xt (y) are conserved quantities under the 
time evolution. 



Lemma 2.1 Assume that {7^, V*} satisfies pd Toda equation H.1Y\1.4\ ■ The 
T\ N , V 1 = T\ N n V t+ 



alld Y\n=l 4 



ntr^ +M - 



Proof. By (& detX t (y) , = - Il„ *2) (Il„ 4 - V> ' (II n v) 

does not change under the time evolution for any y. Then, the set of quantities 
Ut ■= {Un Vn,Un 4, • • • = EL 4 +M ~H satisfies % = U t+ i. From this equation, 
it follows that {UnVlUnll} = {Un ^ + MI„ 4 +M >- The inequality (£3) 
and (|1.4p leads the lemma. ■ 



Let $(x, y) := det(X t (y) — xE) be the characteristic polynomial of X t (y) 
(£? is the unit matrix). The equation 



(2.4) 



defines the affine part C of its completion C. We call this projective curve C 
the spectral curve C of the pd Toda equation. C is a (M + l)-fold ramified 



covering over 



Hi. 



Px ■ C 



(2.5) 



and is also a AT- fold ramified covering p y : C — ► P 1 . It goes without saying that 
C is conserved under the time evolution and is completely determined by the 
initial values {V°, I°,I n ,..., i,f _1 }^ =1 . Note that C may fail to be smooth in 
certain situations. We restrict ourselves to the case where C is smooth. 



Remark 2.1 If M < N, the matrix X t (y) is a N x N matrix and is of the 
form 
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where ok , < i < M , 1 < j < N) are polynomials in {V£, 1^, . . . , 

4 +M_1 }n=i- ^ general, the (i, j)- component of X(y) is the essentially finite 
summation (X(y))ij = X);=-i a 



(j-i+IM+l),, 



y 1 , where a) = Pj> ot 



(M+l) 



= 1 



anda\ P) = (P < -1,P > M+l). 



Remark 2.2 When the greatest common divisor (N,M) ^ 1, i/iis particular 
matrix \2. 6} is not a regular periodic difference operator which is analyzed in 

ED- 

We now analyze the points contained in C \ C. We calculate the polynomial 
expression of the function §(x, y) in x, y and y^ 1 : 

$(x, y) = A (x)y M + A 1 (x)y M - 1 +■■■ + A M {x) + A M+l {x)y- 1 = 0. (2.7) 

For the analysis of the behaviour of C, we need to analyze Aj (x) , (J = 0,1, ... , 
M + 1). The following lemma is established by van Moerbeke and Mumford 
[TTj . They analyzed the determinant det (X(y) — xE) by direct calculation. 

Lemma 2.2 Let m :— (N,M), N — mNi and M — mM\. The polynomial 
Aj(x) is a polynomial of degree kj satisfying 



(0<j<M), k M +i = 0. 



(2.8) 



The equality in \2.8\) holds if and only if the right hand side is an integer. 
Moreover, Aq(x), Amx (x), Aim-l ■ ■ ■ , A m u 1 (x) are expressed as: 



A rMl (^) = (-l) M(JV - M)+r (7) 



.rNi 



We start form the polynomial (|2.7p . Let 7 
obtain the expression 



x N ly -Ah 



(r = 0,1,..., m) 
By lemma 



wo 



y 



-mMi 



V ) = (-l)M(N-M) ( 7 m _ (™) 7 ™-l + (™) 7 ™- 2 + • • • + (-l) m ) 

+ lower order terms when |x| , \y\ — > 00 
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This implies x Nl y 1 = 7 ~ 1 near (x,y) — (00,00). By (N\,Mi) — 1, there 
exists the local coordinate t equipped with the completion C D C such that 

x „ t~ mMl = t~ M , y ~ r mNl = r N , {x, y) ~ (00, 00). 

In particular, there exists only one point P £ C which is expressed as P = 
(00, 00). 

In a similar manner, there exists a point Q €E C which is expressed as 
Q = (oo,0). The local coordinate t near Q satisfies x ~ i -1 , y ~ t . Using 
these fact, the divisors (x), (y) 6 Div(C) are 

(a;) = —MP — Q + (a positive divisor on C), (2-9) 
(y) = -NP + NQ. (2.10) 



Remark 2.3 TTie existence of the unique point P(oo, 00) is a special property 
of X(y) . In fact, there exist m points Pj (x, y) = (00, 00) (J = 1, . . . , ni) on the 
algebraic curve associated to a regular operator matrix \TT$. 

Although the concrete calculations in [TT] should be applied only to the case 
that X{y) is a regular operator, these results are also applicable to our case on 
condition that C is smooth. Precisely, these calculations become true for our 
case by substituting 

P 1 = P 2 = ... = P m (=P). (2.11) 
2.2 The eigenvector mapping 

We now define the isolevel set Tc as the set of matrices X(y) associated with the 
the spectral curve C. The following proposition is fundamental to the algebro- 
geometric method for integrable systems. 

Proposition 2.3 Let C be smooth and X(y) 6 Tq. There is a unique line 
bundle V c C x C N s.t. 

n^ 1 (x,y) = {the eigenspace of X(y) corresponding to the eigenvalue x} C C , 

for all (x, y) 6 C \ {x = 0, 00} where C x C N D V C is called the canonical 
projection. 

By virtue of this proposition, we obtain the map 

Vc- T c -> {U^¥ N - 1 } UCC 
X(y) ~ 

where V v is the dual bundle of V. The section of V v is a component of the 
eigenvector of X(y). By the Grothendieck-Riemann-Roch theorem, it follows 
that 

Im^c C Pic d (C), d = g + N-l, (2.12) 
where g is the genus of C. (See [T2"]). 
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Definition 2.1 For smooth C , the eigenvector mapping associated to the equa- 
tion \2.3\ is the mapping 

<Pc:T c ^ Pic d {C) 
defined as above. We shall call V w the eigenvector bundle. 

The eigenvector mapping is an important tool to analyze the various inte- 
grable systems |13j . The following proposition is essential to our arguments in 
the present paper. 

Proposition 2.4 The eigenvector mapping ipc ■ Tc ^ Pic d {C) is an isomor- 
phism to Inupc ■ 

This proposition is a straightforward result of the following theorem provided 
by van Moerbeke and Mumford. 

Theorem 2.5 There is a one-to-one correspondence between the two sets of 
data: 

(a) a multi diagonal matrix X{y) of the form 12. 6\) such that <&(x, y) = defines 
an affine part of a smooth curve. 

(b) a smooth curve C , two points P, Q, two functions x,y on C and a divisor T> 
which satisfies 

<p c (X(y))=V+(N-l)Q. (2.13) 

„, (JV-l)(M + l)-m+l 

L has genus g = , ana aeg D = g. 

Remark 2.4 Although the equation \2. 13\ does not appear in the van Moerbeke 
and Mumfords paper J 1 If .we easily derive this equation from the relation (p. 107) 

N N-l 

{9k)+V> J2 p i~ X)Q<> k = l,2,...,N-l, (2.14) 

i=k-\-\ i—k 

where (g\, . . . , <?at— 1> 1) T is a section of the eigenvector bundle V y . In fact, 
prm > yields 

(9k)ac < T> + (N - k)Q, (2.15) 

which implies d — deg (hmpc) < degP + iV — 1 = g + TV — 1. Because of the 
equality $2.12\ , (31)00 must satisfy 

(ffi)oc =V+(N-l)Q. (2.16) 
(fQ5|) and (fO^|) imply ifQ^l . 

Definition 2.2 A finite component (p{ n (X(y)) of the eigenvector mapping is a 
positive divisor T> of degree g which appears in \2.13\i . 
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2.3 The action of the evolutions on the eigenvector bundle 

In this section, we represent the two actions — index evolution (n i— ► n+1) and 
time evolution (t >— » t + 1) — on the eigenvector bundle. 

Proposition 2.6 Let D n be the divisor D n = P — Q. Then the following dia- 
gram is commutative. 



n\— m+1 -1 



Pic d {C) 
Pic d {C) 



Proof. Let us denote a : n t— ► n + 1. A straightforward calculation leads to 



X(y) 



( V1 \ 




( V1 \ 




= X 




Vn-i 




UJV-1 


\ v N J 







/ V 1 vn\ ( V 1 v N \ 



« ^(A^,)) 



\ WJV-1 / 



(2.17) 

It is enough to prove (y 1 g N 1 _ 1 ) c-(gi)o ~ ~D n . By we have (y5w-i)+ 

D > (A — 1)Q — (A — 1)P. An argument similar to remark l2~^l allows to conclude 
(j/ _1 .9tv-i)oo ~ (yffjv-i)oc = f + (A — 1)P, which completes the proof. ■ 

In order to determine the action of the time evolution on the eigenvector 
bundle, we introduce the concepts of Bloch solution and transposed operator. 

We identify the eigenvectors of X(y) with the Bloch solutions with multiplic- 
ity y of the periodic infinite matrix 



X 



\ 



(i) 



'at 



(i N a\ 




„(2) 



01 «2 1} 





AM) 
*M-l 



J 2 ) 



1 

v (M) 
M 





1 





1 



V 



which are the infinite vectors v — (■ ■ ■ , f n -i, w„, v n +i, • • • ) T such that 

Xv = xv and v n+ N—yv n . (2-18) 
The first equation of (]2 . 18[) is equivalent to 



Pn-lVn-l + a^'Vn + a^]_ x V n+ x H h a\*+M-iVn+M-l + v n+M = xv» (2.19) 

for n £ Z (o^+at = a ™^> /Wa 7 — Ai)- Because the l.h.s. of (|2.19|) is a lin- 
ear combination of u n -i, v n , ■ ■ ■ ,v n+ M-i, the Bloch solution associated with 



(M) 
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(x, y) € C are also be written as a linear combination of v^ 2 \ . . . , v^ M+1 \ 
where 



V' - ( 10 v ( ' v ( ' V 

< — y. . . , i, u, . . . , u, u M+2 , u M+3' " ' ) 

(2) _ / o i o (2) (2) \T 



w M+2> 

SjCW+l) _ f Q Q ! (M+l) , T 

More precisely, let ip(x,y) be the Bloch solution: 

i/j{x,y) = aiSr^ + 02«r 2 ^ H h a M +ii> (M+1) 



(2.20) 



where = a,i(x,y) and tr*' = t/ l ^(a;). 

Recalling (|2.3[) and proposition ^. 31 the eigenvector t>t(a;, y) at time t satisfies 



v t+1 (x,y) = R t (y) ■ v t (x,y). 
Equivalently, the Bloch solution ip t (x,y) satisfies 

ij t+1 (x,y) = Rf^(x,y), 



(2.21) 



(2.22) 



where 



Rt 



N 1 



It 1 



ft 



Tl 



(j2T2"0"j) and (|2~22]) yield 



,t+i 



a\(R t v {1) ) j + --- + a t M+1 (R t v 



(M+l), 



j = 1,2,..., M+l, (2.23) 



where (v)j is a jth-component of the vector v. Equation (12.23)) is equivalent to 



/ t + l t+l t + l t+l "\ T _ TT t t t t t \ T 

\ a l , a 2 :---: a M '"M+lJ — a t ' \ a l> "2> • • • ) "M> a M+l) J 



where 



/ 1 



H t = 



1 2 



(1) (2) 
A/+2 ^M+2 



Tl 

1 M 



(2.24) 



. . . V (M) p , (M+l) / 
u Af+2 Z M+1 ^ ^m+2 / 



Using this equation, we obtain 

[v{ + \v t +\ w&Si&kf = H t ■ (v\,vi v* M , vi I+1 ) T , (2.25) 
where ^{x^y) = (. . . , v\ ,v\, . . . , v* M+l , . . . ) T . 

Lemma 2.7 For fixed generic x, the M + 1 Block solutions associated with x 
are linearly independent. 

Proof. For generic x, associated multiplicities yj (j = 1, 2, . . . , M + 1) of the 
Bloch solutions are all distinct. ■ 

Lemma 2.8 detif t = {-l) M+1 I\x. 

This lemma is proved by an elementary calculation, which we shall give in the 
appendix. 

Let T be the invertible sheaf associated with the eigenvector bundle V v . 
Let us consider the direct image (p x )*J-, where p x is a projection of C defined 
by (|2.5[) . The sheaf (px)*^ is a locally free sheaf of rank M + 1. By lemma [2~7l 
the fiber of this direct sheaf is C M+1 . 

Equation ()2.25|) is regarded as the time action to the space C M+1 . Re- 
call that the components of the Bloch solutions of X are the section of the 
sheaf T. For x G P 1 , we denote the fe-th component of the infinite vec- 
tor ^} l {x,yf) by Vj k (x), and ^(x^j) by ipj(x). By (|2.20p . the finite vectors 

i>) ■■= v] a , v t jM+1 ) T have the property 

Wjljil 1 are linearly independent {ipjjf^ 1 are linearly independent. 

(2.26) 

Equation (|2.25[) implies that 

• ■ • , fi A V + i) =H t ft M+1 ). (2.27) 

On the other hand, lemma 12. 71 12.81 and (|2.26[) imply that there exists at least 
one vector -ip* which satisfies 

mult ;r -0* +:L > mult^j, (2.28) 

where mult^(zi, Z2, ■ ■ ■ , zn) T — niin [ mult x zi, mult;^, . . . , xmx\t x ZN ] and 
multzZ is multiplicity of x in z. 

The same discussion can be repeated for the projection 

Py :C^¥\ 

which is an ./V-fold ramified covering over P 1 . 

The following two facts are then obvious to prove. 

Lemma 2.9 For fixed generic y, the N eigenvectors of Xt(y) are linearly in- 
dependent. 
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Lemma 2.10 det Rt{y) is a polynomial of degree one in y. 

Let us consider the direct image {p v )*J-. From lemma HOfl we then find that 
the fiber of this direct image is C N . 

Equation (|2.2ip can be regarded as the time action to the space C^. For 
fixed y, we denote p y 1 {x) — {{x\, y), . . . , (xn, y)}. Then we obtain that 

3j s.t. mult(y- yt) v t +i(xj,y) > mu\t (y _ yt )Vt(xj ,y), (2.29) 

where det R t (y t ) = y t = Xljl). 

2.4 The transposed operator 

In this subsection, we introduce the transposed operator, and give the proof of 
proposition 12.161 

The section of the eigenvector bundle V v can be described as a rational 
function of x and y. We denote the set of these sections by F(V V ). Let Aij := 
(-1)*+*' x (i, j)-th minor of X(y) -xE. We have g k = ^{k = 1, 2, . . . , N- 1), 

where (gi, . . . ,gN-i, 1) T 6 T(V V ). We are interested in the divisor (51)00 and 
hence, it is important to explore the common zeros of Ajy,i and Ajv,jv- 
Equation (|2.3p is equivalent to 

X t (y) T R t (y) T = R t (y) T X t+1 (y) T . (2.30) 

Let us define s := —t, and A* := JA T J for a regular matrix A where 

/0 1\ 

J = 1 • 
\1 oj 

Then (|2~30|) becomes 

X* +1 R* +1 = R* +1 X*. (2.31) 

Note that the matrix X*(y) is also of the form (|2.6p . We call this new matrix 
X* the transposed operator of X. A careful analysis of the eigenvector bundle 
of the transposed operator gives more information on the original matrix X. 

Lemma 2.11 (van Moerbeke, Mumford) (1) There exist positive regular 
divisors T>\ , T>i , and P3 such that 

(A N ,i/A N , N ) = V 1 + (N- 1)P - <p {n (X) -(N- 1)Q, (2.32) 

and 

{A ltN /A N , N ) =T> 2 + (N - 1)Q - V 3 -(N- l)P. (2.33) 

Moreover, degPj = g (i = 1,2,3). 
(2)The divisor (Ajv,jv) satisfies 

(A N>N ) = cp in (X) +V 3 - (NM -M-n+ 1)P - NQ. (2.34) 
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Remark 2.5 Lemma [2.11\ (l) and Remark ] 2. 16\ imvlv that dee 
1, and that T>2 and P3 do not have common points. 

Lemma 2.12 <pfo(X*) = ip in {X T ) = V 2 . 



= g + N- 



Proof. Note that X* and X T give the same spectral curve C . By definition, it 
follows that: X T (Jv) = x(Jv) X*v = xv, which implies the first equality 
of the lemma. 

The second equality is obtained from (|2 . 13[) . (|2.33[) and the fact that the (N, 1)- 
th minor of X is the (1, N)-th minor of X T . ■ 

Lemma l2~12l and (|2~33l) yield 

V 3 ~ tfaiX*) + (N — l)Q -{N- 1)P ~ <ff a (X*) -Q + P. 

Using proposition 12.61 we obtain D3 ~ (pf n (a(X*)) = V5f n ((cr _1 X)*) . In fact, by 
virtue of the Riemann-Roch theorem, we obtain the following stronger result: 

2?3 = ^((a^X)*) 

because T>z is a regular divisor of degree g. 

Proposition 2.13 Ajv.jv has 2g zeros in C\{x = 00}. And the corresponding 
divisor is 

i Pln (X) + i Pfn ((a- 1 X)*). 

As a corollary of this proposition, the zeros of A^i, A^jv, Ajv,i are also deter- 
mined: 

Corollary 2.14 A1.1, Ai.jv, Ajv,i have 2g zeros in C\{x — 00} respectively. 
The corresponding divisors are: 



A M 

Ai.at 
A ATI 



WnicT^X) + <p in (X*) 
<Pfn(X) + <pUX*) 

ifUa^X) + ipUia^X)*). 



(2.35) 



Proof. The divisor of A1.1 is obtained from Ajv.jv by the change of the index: 
n — > n — 1. The remaining divisors are determined by the identity Ai^Ajv.jv = 
Ai,jvAjv,i and (|2T6]l . ■ 

Example (N = 4, M = 2) 
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The correspondence X <-* (er X)* is equivalent to the correspondence: 

(1) (1) (2) (2) a a 

Our aim is to analyze the distribution of common zeros of Ajv.i and Ajv,iv- 
This set of zeros has the following characteristic property. 

Lemma 2.15 

{ The common zeros of {Ajv.i, A^jv}} 
= {The common zeros of {Ajv,fc}fe=i,2,...,jv}- 

Proof. Consider the following open subset on C: C = C\{x = oo}. By (|2.15[) 
and (|2.16p , for k > 2, we obtain (gk)oo\c < {gi)oc\jj- This is equivalent to 
(Ajy.fc/A^iv)^ (Ajv.i/Ajv.Ar)^ , which proves the lemma. ■ 
Using the preceding calculations, we come to the linearization result. 

Proposition 2.16 Let D(j\ (j = 1,2, ...,M) 6e i/ie divisors Dtj\ = Aj — Q, 
where Aj = (0,yj) and yj is the complex number which satisfies y,j = Yi n ^n- V 
t = j (mod M), the following diagram is commutative. 

T c -> Pic d {C) 

ti->t+l I I +D {j) 

T c -> Pic d (C) 

Proof. Recall that ]\ n T n = n„ J n +M b y lemma O Let t = j (mod M). 
Proposition 12.131 and corollary 12.141 imply 

{Common zeros of Ajv,i and Ajv.jv} = ffn((o'^ 1 X)*). 

Using this fact, it becomes clear that the equations (|2.28[) and (]2 . 29[) give the 
time evolution action of the transposed operator X* . In fact, by virtue of lemma 
12.151 (|2.28[) requires the existence of a point (x, y) = (0, y') which becomes the 
common zero of Ajv,i and Ajv,jv when t 1— * t + 1. Then (|2.29|) also requires the 
existence of a point (x,y) = (x' ,yj) which becomes the common zero of Ajy,i 
and Ajv.jv- Using the same argument about poles, it follows that the divisor 
V' [rj) = Q-Aj obeys the relation: ip {n (X^ +1 ) = ip in (Xt)+V\j). Recalling (f2T3l]) . 
the two time evolutions defined by X t and X* are opposite: s = — t. Hence the 
divisor — —T^'(j\ = Aj — Q gives the time evolution of original system. ■ 

3 Theta function solutions 

3.1 The distribution of the points ipfn(X) 

Let <£>f„ (X) = Pi H \-P g , Pj = {xj ,yj) € C. We are interested in the numbers 

Xj E C, (j = 1,2, ... ,g). Now, consider the resultant of the polynomials [14] 
$>(x,y) = ydet (X(y)-xE) and Ajv.iv as polynomials in y. Let us denote this 
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resultant by Rest y (<£>, Ajv,jv) —■ R(x), which is a polynomial in x. More pre- 
cisely, R(x) is an element of C[a^ 1 ' ) , . . . , a$ ; . . . ; ok^ , . . . , a^ 1 ^ ; /3i, . . . , @n] [x]. 
From proposition ^. 131 it follows that deg^, R(x) = 2g. 

We also consider the resultant S(x) := Rest a ($, A^jy). Due to proposition 
12.131 and corollary 12.141 the common divisor of R(x) and S(x) is an element 

of C(a^ , . . . , , . . . ; , . . . , /3i, . . . , /3jv)[ie], the degree of which (as a 
polynomial in x) is equal to g. Multiplying divisors (if needed), we obtain the 
monic polynomial 

T(x) £ C(a^\ . . .,a$; . . . ;a^ M) , . . . , a^f ; /3i , . . .,(3 N )[x], 

the roots of which are the common roots of R(x) and S(x), i.e., the common 
roots of A^jv and Ajv,jv- Recalling that the set of the common zeros of Ajv.jv 
and Aijv contained in C \ {x — oo} is {Pi, . . . , P g }, we conclude: 

deg, T(x) = 5 , T(^) = 0, i=l,2,...,g. (3.1) 
3.2 Theta function solutions 

For a complex curve (or Riemann surface) C of genus g, one usually considers 
a canonical basis of H\(C,Z). We denote the canonical basis by a%, . . . ,a g ; 
bi, . . . , b g £ H\{C, Z). Let uj\, . . . , u> g be the holomorphic differential of C which 
satisfies J a .^i — 5j,i- A period matrix of the Riemann surface C is a g X g 
matrix £> = (Jl.Wj). Let 9(z,B) be the t/ieta function; C 9 — > C, and A : 

Pic 9 J(C) (: = C S /(Z» + BZ 9 )) the 4&eZian mapping. The following theorem 
is a classical and fundamental result. 

Theorem 3.1 (Riemann) Let C be a Riemann surface of genus g, and let 
T> = Pi + • ■ ■ + P g be a regular positive divisor. Then the function 

F(p) = 8(A(p) - A(V) — K,B), P £C 

has exactly g zeros p = Pi , . . . , P g on C , where K is the Riemann constant of 
C. 

To obtain the solution to the pd Toda equation, we consider the following inte- 
gral: 

2m J 9C o X ^ F(p) ( ' 2m j dC o dZ ( p ^j ' ^ 3 ' 2 ' ) 

Here the integral path dC° goes along the edge of the simply connected domain 
C° obtained from the Riemann surface by cutting it along ai, . . . , a g ; bi, . . . , b g . 
The integral I can be rewritten as 

i - & ELi {k + + k + V) dz( P ) 
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Recalling the classical fact F(p + a&) = F{p) 1 F{p + bk) = cxp(— 2Tri(A(p) — 
A(V) — K)k)F(p) and d(A(p))k — ^k(p), the integral / is transformed to 

I = ELlL k X (P)^( P ). (3-3) 

On the other hand, by the residue theorem, the integral / also has the expression: 
1 = Ef=i x(Pi) + Resp(dZ) + Res Q {dZ). (3.4) 
Let tp and £q be local coordinates around P and Q respectively. These satisfy 

x ~ l/(tp) M , (neighbor of P), x ^ l/tg, (neighbor of Q). 
In a neighborhood of P, one has: 



(t P ) M dtp F (t P ) M ^ v & '\ dtp ) 



(3.5) 



To calculate the residue of the differential (|3.5[) . we explore the behaviour of 

d< d±~p l — w i around points P and Q. Let q :— Resp(cj//(ip) M ). Then we 
obtain the expression: 

Resp(dZ) = £f =1 log(P(P))). (3.6) 

In the similar manner, we also conclude 

Res Q (dZ) = Ef =1 c{(ft log(P(Q))), (3.7) 

where c[ := ResQ (w; /ig ) . 

By ([331) and (33|), we obtain 



X>(^) = E /^(p)^b)-E c '^ lo g( F ( p )))-E^^ lo g™)))- (3-8) 

Using this equation, we obtain the following theorem which is a generalization 
of the preceding result concerning the theta function solution to the pd Toda 
(M = 1) equation 10J and which is the main theorem in the present paper: 

Theorem 3.2 Let T(x) be the monic polynomial of degree g obtained by <\3. 11 : 
T(x) = x s - aix 9 ' 1 + ■■■ + (-1)%, 

with a%,...,a g £ C(a^, . . . ,or^,- ■ ■ ;a[ M \ . . . , ; /3i , . . . ,/3jv). Then we 
find 

a 1 = J2Li I ai x (p) "lip) -YlUiCi {d^og6(nk + u{t) + c a , B)} 

~ Ef=i 4 {di log0((n + l)fc + u{t) + c , B)}, (3.9) 

where c ; = Res P {ui/ {t P ) M ), c[ = Res Q {uJi/t Q ), k = A(P - Q), c = A(Q - 
V Q + e), andv{pM + q) = pA{A x + - ■ ■ + A M ) + A(A 1 + - ■ ■+A q )-(pM +q)A(Q) 
(1 < q < M). Here is the theta divisor: A(Q) = — K . The divisor T>q is the 
initial value tph\{Xt=o{y)) = T>q. 
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A Proof of lemma 12.81 

In this section, we give the proof of lemma 12751 Let 

U k := R t+k -i(y)Rt+k-2{y) ■ ■ - Mv), l<k<M. 

Define the row vector «( fc ) = (u^tifj, . . . ,u^ N ), where Uf. = (u^j ), and the 
following homomorphism of rings: 

a : Z [{Vl 7*}„ ez ] - Z [{V; 4 , /*} neZ ] ; V£ -» 7* -> 

We rewrite the j'-th component of as Uj for short. By definition of Rt(y), 
we obtain 

+1) = a(uf\) + l\+ k uf\ («W = 0). (A.l) 
On the other hand, the 2nd row of the matrix X t (y), which is of the form 

(/?i,4 1 \a£ 2 \... ) a£ +1) l,0,...,0), 

satisfies /?i = V/u^, a< j+i = v{ u ^^) + Vi u j+i- Figure Q] displays the al- 
gorithm which we shall use to obtain the row vector expressed by (|A.1[) . 
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/X /\ /\ /x 

Figure 1: The j + 1-th row of this diagram displays the first j + 1 non-zero 
components of the row vector m®. 

Let us introduce the signs S and \ to describe the terms which appear 
in the figure [TJ Let us define the set of arrows Ar := and the set of 

sequences Ar r := {(ai, a 2 , . . . ,a r ) \ Oj 6 ^4r}, r£R We define the map of sets 
{■} : U r ArT Z[{/*+ / } 1 <„< w , 1 < / < A/ ] as follows: 
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For r = 0, define {0} := 1. If {a\,a%, ■ ■ ■ , a r } S Z^i^" 1- '},^] is given, we define 
{oi, a 2 , . . . , a r , a r+ i} by 



,/}:=^+ r { ai ,a 2 , 



4, {ai,...,a r ,\} :=cr({ai, 



-}) 



inductively. For example, {/} = 7*, {//} = {/\} = /*, {\/} 

J* +1 . By definition, the j-th component of satisfies 



,(*0 



{ai,...,a fc }. 

It\=j-l> $S=k-j+i 



(A.2) 



Lemma A.l {/, a 2 , . . . , a fc } = ijf +1 {\, a 2 . . . , a fc }, 
where I = jj{i/ie arrow \ included in {a 2 , . . . , afc}}. 

Proof. We prove the equation by induction respected to fc. If k = 1, the 
equation is equivalent to {>/} = -f*{\} which is true by definition. Let fc > 2. 
If afc = /, we obtain l.h.s = I^~ x {/, a 2 , . . . , a k ^}, r.h.s = /^^{N 
,a 2 , . . . ,afc_!}. By assumption of induction, it follows that l.h.s = r.h.s. We 
also prove the equation in the similar manner if —\. H 



Lemma A.2 u[ M) + Y,f=i • • -I) = 0. 

Proof. 

4 m) +e;i 1 (-i)^i--^ 



= {//---/}+E(- 1 ) 3 



3=1 
M-1 



tt\=3, l^=M-i 



= {//••■/}+ E E •!/;•••/;/ 

3 = 1 S\=3+l,l)^=M-j-l 

M 

+ E(" 1 )' Ei^'*}^-" 7 3 (•••LemmaEI 

= {//■■■ /}-h{\/ ■■■ /} 

= ■ 



t 

3+1 



ft . . . J* 
J l i 3 



Lemma A.3 E^i (-l)^^)^ •••/] = {-1) M I{ ■ ■ ■ I M+V 
Proof. We start from lemma IA~2l 

' AM) 



=^r ) )+E^ 1 (-i)M4+i ) )^--- j i + i- 
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Multiplying I\, we obtain 

o = E i =l 1 (-iM«?°)^-^ 
= E-U (-i)v(^ M) )/l • • -3 + (-i) M+1 /f • • • 

which complete the proof. ■ 
proof of lemma 12.81 

Calculate det H t by the definition of H t (|2.24p . The components 4+2 7 (j = 
1, 2, . . . , M + 1) are rewritten as 

4+2 = "A, 4+2 = X ~ > 4+2 = -^"^ > CJ = 3, 4, . . . , M + 1). 

by virtue of (|2.19p . The expansion of the determinant with respect to the 
(M + l)-st row yields 

det H t = (-l) M {-/3i + (a^ - x)I{ - af + • ■ ■ 

=(-i) M {-th-xn - Ejli • • • 3 + i-v M n ■ ■ ■ ^Wi} 
= (-i)« {-^4 M) - ^ - s£i (-lyMu? ) + *?«S5m • • • 7 i 

+(-i) M ii---4ii f+1 } 

= {-l) M+1 l\x. {•; Lemma [O] and \EM ■ 
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